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Abstract. We provide two methods for constructing smooth bump functions 
and for smoothly cutting off smooth functions on fractals, one using a prob- 
abilistic approach and sub-Gaussian estimates for the heat operator, and the 
other using the analytic theory for p.c.f. fractals and a fixed point argument. 
The heat semigroup (probabilistic) method is applicable to a more general class 
of metric measure spaces with Laplacian, including certain infinitely ramified 
fractals, however the cut off technique involves some loss in smoothness. From 
the analytic approach we establish a Borel theorem for p.c.f. fractals, showing 
that to any prescribed jet at a junction point there is a smooth function with 
that jet. As a consequence we prove that on p.c.f. fractals smooth functions 
may be cut off with no loss of smoothness, and thus can be smoothly decom- 
posed subordinate to an open cover. The latter result provides a replacement 
for classical partition of unity arguments in the p.c.f. fractal setting. 



1. Introduction 

Recent years have seen considerable developments in the theory of analysis on 
certain fractal sets from both probabilistic and analytic viewpoints [TJ Q3J [55] . In 
this theory, either a Dirichlet energy form or a diffusion on the fractal is used to 
construct a weak Laplacian with respect to an appropriate measure, and thereby 
to define smooth functions. As a result the Laplacian eigenfunctions are well un- 
derstood, but we have little knowledge of other basic smooth functions except in 
the case where the fractal is the Sierpinski Gasket [TH [5J HHj ■ At the same time 
the existence of a rich collection of smooth functions is crucial to several aspects 
of classical analysis, where tools like smooth partitions of unity, test functions and 
mollifications are frequently used. In this work we give two proofs of the existence 
of smooth bump functions on fractals, one taking the probabilistic and the other 
the analytic approach. The probabilistic result (Theorem [^H]) is valid provided the 
fractal supports a heat operator with continuous kernel and sub-Gaussian bounds, 
as is known to be the case for many interesting examples [TJ [3] that include non- 
post-critically finite (non-p.c.f.) fractals such as certain Sierpinski carpets. In this 
setting it can also be used to cut off functions of a certain smoothness class, with 
some loss of smoothness (Theorem l2.7l) . By contrast the analytic method (Theorem 
3.8|) is applicable to the smaller class of self-similar p.c.f. fractals with a regular 
harmonic structure and Dirichlet energy in the sense of Kigami [14] . 

For p.c.f. fractals we use our result on the existence of bump functions to prove a 
Borel-type theorem, showing that there are compactly supported smooth functions 
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with prescribed jet at a junction point (Theorem 14. 3|) . This gives a very general 
answer to a question raised in [19( [5] , and previously solved only for the Sierpinski 
Gasket [3D] . We remark, however, that even in this special case the results of [2U] 
neither contain nor are contained in the theorem proven here, as the functions in 
[20] do not have compact support, while those here do not deal with the tangential 
derivatives at a junction point. 

Finally we apply our Borel theorem to the problem of partitioning smooth func- 
tions. Multiplication does not generally preserve smoothness in the fractal setting 
[3], so the usual partition of unity method is not available. As a substitute for this 
classical tool we show that a smooth function can be partitioned into smooth pieces 
with supports subordinate to a given open cover (Thcorcm l5.ip . This result plays a 
major role in a forthcoming paper on a theory of distributions for p.c.f. self-similar 
fractals [22] . 

Setting. Let X be a self-similar subset of M. d (or more generally any complete 
metric space) in the sense that there is a finite collection of contractive similarities 
{Fj}jLi °f the space and X is the unique compact set satisfying X — UjL 1 Fj(X). 
The sets Fj(X) are the 1-cells, and for a word w — (w\, u>2, ■ ■ ■ , w rn ) of length m we 
define F w = F Wl o • ■ ■ o F Wm and call F W (X) an m-cell. If w is an infinite word then 
we define [w] m to be its length m truncation and let F W (X) = n m Ft w ] m (X), which 
is clearly a point in X. The map from infinite words to X is surjective but not 
injective, and the points of non-injectivity play an important role in understanding 
the connectivity properties of the fractal (see Section 1.6 of [13] )■ In particular 
there are critical points of the cover by the Fj, namely those x and y for which 
there are j ^ k in {1, . . . , N} such that Fj(x) = Fk(y) (so Fj(x) is a critical value). 
We call an infinite word w critical if F W (X) is a critical value, and then call w 
post-critical if there is j e {1, ... , N} such that jw is critical. The boundary dX of 
X consists of all points F^,(X) for which w is post-critical. In the case that the set 
of post-critical words is finite the fractal is called post-critically finite (p.c.f.) and 
we also use the notations Vq = dX and V m — U w F w (V~o), where the union is over 
all words of length m. The points in (U rn V m ) \ Vq are called junction points. We 
shall always assume that Vq contains at least two elements. 

We suppose that X comes equipped with a self-similar probability measure 
fi, meaning that there are /ii,...,/ijv such that the cell corresponding to w = 
(wi, . . . , w m ) has measure fi(F w (X)) — Yij=i f-wj ■ In order to do analysis on X we 
assume that X admits a Dirichlet form £, so £ is a closed quadratic form on L 2 (/i) 
with the (Markov) property that if u £ dom(£) then so is u — uxo<u<i +x«>i an d 
£(u,u) < £{u,u), where xa is the characteristic function of A. We will work only 
with self-similar Dirichlet forms, having the property that 

(1.1) £{u,v)= r- 1 £{uoF w ,voF w ) 

m-words w 

where the factors rj are called resistance renormalization factors and as usual r w — 
i~wl ' ' ' r w m ■ For convenience we restrict to the case of regular harmonic structures, 
in which < r j < 1 for all j. In addition we assume £ has the property that 
C(X)ndom(£) is dense both in dom(£) with £-norm and in the space of continuous 
functions C(X) with supremum norm. We often refer to £ as the energy. If X is 
a nested fractal in the sense of Lindstr0m [17] then such a Dirichlet form may be 
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constructed using a diffusion or a harmonic structure [El [23] . Other approaches 
may be found in [21] H3 UM 03] • 

Using the energy and measure we produce a weak Laplacian by defining / = 
Au if £{u, v) = —Jfvd/j, for all v G dom(£) that vanish on dX. Then —A 
is a non- negative self-adjoint operator on L 2 (/j,). When Au G C(X) we write 
u G dom(A); this notation is continued inductively to define dom(A fe ) for each k 
and then dom(A°°) = rifedom(A fe ). We say / is smooth if / G dom(A°°). On 
a p.c.f. fractal the weak Laplacian admits an additional pointwise description in 
which £ is a renormalized limit of energies £ m corresponding to the finite graph 
approximations V m and the Laplacian A is a renormalized limit of the associated 
graph Laplacians A m . Details are in Section 3.7 of [14] . 

By standard results, existence of the Dirichlet form £ implies existence of a 
strongly continuous semigroup {Pt} with generator — A. Conversely if there is such 
a semigroup and it is self-adjoint then there is a corresponding Dirichlet form, so 
we could equally well begin with {Pt} and construct £ (see OH]). The Markov 
property of £ ensures that if < u < 1 fi-a..e. then also < Ptu < 1 /i-a.e. 

If X is p.c.f. then there is a definition of boundary normal derivatives of a 
function in dom(A) and a Gauss-Green formula relating these to the integral of 
the Laplacian on X . The usual definition uses resistance-renormalized limits of the 
terms of the graph Laplacian that exist at the boundary point. If qi is the boundary 
point of X that is fixed by Fi and ri is the resistance factor corresponding to F- L 
we may define a normal derivative d n at qi and have a Gauss-Green formula as in 
Section 3.7 of [H] by 

(1.2) d n u( qi ) = - lim r^ m A m . q% u{ qi ) 

in—*oo 

(1.3) (v(q)d n u(q) - u(q)d n v(q)) = / (uAm - uAi))d/j 
qeax J x 

where in (|1.2[) the quantity A mtqi u(qi) is the graph Laplacian at qi. (We distinguish 
this from the graph Laplacian at interior points because it may involve a sum over 
fewer edges - for example in the case of the Sierpinski gasket the graph Laplacian 
at interior vertices has four terms, but at boundary points there are only two.) 

Normal derivatives may also be localized to cells, so that d nj F ul (X)u{F w {qi)) is 
given by the limit in (|1.2|) but with A m ^ w u(F w (qi)) denoting the graph Laplacian 
restricted to to edges that lie inside F W (X). It is then easy to see that if Au 
exists and is continuous on each of finitely many cells that meet at F w (qi), then 
it is continuous on their union if and only if the following conditions hold: u is 
continuous, Au has a unique limit at F w (qi), and the normal derivatives at F w (qi) 
sum to zero. We call these the matching conditions for the Laplacian. 

We shall have need of two other pieces of information about a p.c.f. fractal with 
regular harmonic structure. The first is that there is a Green's function g(x,y) > 
that is continuous on X x X and has self-similar structure related to the discrete 
Greens function ^(x,y) on (Vi \ Vq) . According to Section 3.5 of [14] 

m— 1 

(1-4) 9(x,y)= lim ]T r w ^(F~ 1 (x),F- 1 (y)) 

k=0 weW k 

where Wk is the collection of words of length k but the sum is only over those w 
such that F~ 1 (x) and F^ 1 (y) make sense. Integration against (the negative of) 
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this Green's function gives the Green's operator Gf(x) = — J g(x,y)f(y)dfi(y) and 
solves AGf — f with Dirichlet boundary conditions. In particular we will make use 
of pointwise estimates of g(x,y) that follow from (|1.4|) . The second thing we need 
to know is an estimate on the oscillation of a harmonic function on a cell F W (X), 
details of which are in Section 3.2 and Appendix A of |14j . A harmonic function h(x) 
is determined by its values on the boundary Vo, and its values may be computed 
using harmonic extension matrices Ai via h\F w (v ) = ^w m ' ' ' -^wi ' ^ V6 ■ The A4 are 
positive definite, have eigenvalue 1 on the constant function and second eigenvalue 
at most rj. It follows immediately that the oscillation of Aj/i|Vb is at most rj when 
the oscillation of h\v is bounded by 1, and similarly that the oscillation of h\jr m (x) 
is at most r w . 

More details about analysis on self-similar p.c.f. fractals may be found in [M] . 
or [26] in the special case of the Sierpinski Gasket. The lecture notes of Barlow p] 
cover the probabilistic approach that begins with a diffusion semigroup; non-p.c.f. 
examples include Sierpinski carpets 2 , 3J . Some of the general theory connecting 
Dirichlet forms, heat semigroups and spectral theory of the Laplacian is covered in 



Smooth bump functions. In classical analysis on Euclidean spaces the usual 
bump functions to consider are of the form u G C°° with support in a specified 
open set fi, bounds < u < 1 and the property u = 1 on a specified compact 
K C O. In the fractal setting described above it is not usually the case that a 
product of smooth functions is itself smooth (see [4] , or Section [5] below) , so there 
is less practical benefit to asking that our bump functions be identically 1 on X 
and we will not always do so. Nor is it always essential that < u < 1, though this 
is sometimes useful. For this reason we will use the term smooth bump function 
to mean a function u G dom(A°°) with support in a specified set Q and a bound 
\u — 1 < e on a specified compact K C f2. 

Suppose A is a p.c.f. self-similar fractal and we have a function u G dom(A°°) 
with \u — 1| < e on K C A and such that A k u and d n A k u vanish on Vq for all k. 
Then for any word w we see that 



is a smooth bump function with support in F W (X) by the matching conditions for 
the Laplacian. For this reason we also use the term smooth bump function to refer 
to u G dom(A°°) on X with \u— 1| < e on K and which vanishes to infinite order 
at all points q G Vq, by which last phrase we mean A k u(q) — d n A k u(q) = for all 



In this section (A, dist) is a metric space with Borel measure /1 and a self-adjoint 
Neumann Laplacian A. Our main result here is a procedure for using the heat flow 
on X to cut-off a smooth function with only a small loss in smoothness. For this 
purpose we require two assumptions on A. 

The first assumption is that A has a positive spectral gap, in the sense that there 
is A > such that the spectrum of A is contained in {0} U [A, 00). This implies the 
estimate \\P t — IW2.2 < min{At, 2}, where || ■ 1 1 2,2 refers to the operator norm on L 2 , 



mm- 




on the cell F W (X) 
elsewhere 



k. 



2. A SMOOTH BUMP FROM THE HEAT OPERATOR 
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Pt is the heat operator at time t and / is the identity. For future use we also note 
that the spectral representation Pt = / °° e~ xt dE&(x) implies there is an estimate 

(2.1) l|A fc P 4 |k 2 <c k r k . 

For the second assumption, define 



(D{t,d)) 2 = supl / \P t f\ 2 : \\f\\l — 1 and Lcl with dist(Sppt(/), L) > d 




so that for any set L and any f £ L 2 with dist(Sppt(/), L) > d we have 
(2-2) {l L \ Ptf \ 2 ) ' ^MH/Ha- 

Our second assumption is that there is eo so that for any < e < eo there is a 
decreasing sequence tj with all tj < j and £V tj = T < oo, and such that for every 
non-negative integer fc, 

oo 

(2.3) t 'J+i D ( t J^^ 3 ) = Ck < oo. 

3=1 

It is worth noting that if this condition is satisfied then we may make T > as 
small as we like, because D(t, di) < D(t, c?2) if > di and therefore 

oo oo 

J2t7+ m+ iD(tj+m,e2-i) < ^t7 + fc m+1 P(< J+m ,e2-^+ m )) 

J=l 3=1 

OO 

= E tf k +iD(t ri e2-i') <C k 

j'=m 

so we may simply take m so that J^T tj+m ^ s as small as desired and use the 
sequence tj = tj +m . At the same time, this allows us to reduce a finite number of 
the Ck to be as small as we desire. We suppose Co < \. 

The second assumption is given in this somewhat artificial form so as to empha- 
size the estimates needed in the proof. It is a non-trivial assumption, but is known 
to be true in many examples. For instance, if Ptf is given by integration against a 
kernel p(t, x, y) that satisfies a sub-Gaussian upper bound 

, N 7i / /dist(a;,i;) /5 \ 1 /(/5-i)\ 

(2.4) p( ^y)<Jl_exp^- 72 ( yi^J J 

then (|2.3p may readily verified for sequences tj that decrease sufficiently rapidly, for 
example tj = j~i or tj — e~i 2 . Heat kernel bounds like (|2.4[) have been the subject 
of a great deal of research, not only on fractals but also on manifolds and graphs 
(see [SJ [Ml EH E] and the references therein) . Here we satisfy ourselves with noting 
that they are known for many interesting examples, including various p.c.f. fractals 
[TBI [3 [23l [12] and certain highly symmetric generalized Sierpinski carpets (which 
are not p.c.f.) [2, 3J. They are also valid on products of some fractals 25J. Later we 
shall see that requiring p(t, x, y) be continuous and satisfy (|2.4[) is an appropriate 
choice when using our cutoff procedure to produce a smooth function. 
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Definition 2.1. The space of heat-smoothed functions A(X) on X is 

A(X)= \Jp t (L 2 (X)) 

T>0 

Functions in A are characterized by the property that when expanded with respect 
to an orthonormal basis of Laplacian eigenfunctions, the coefficients decay expo- 
nentially with respect to the sequence of eigenvalues. In particular A(X) contains 
the eigenfunctions and is thus dense in L 2 {X). 

Theorem 2.2. Let 4> e A(X) be cj> = P T f for some f € L 2 (X) with ||/|| 2 < 1, and 
T > 0. Let K be a fixed compact set. IfO < e < eo then there is a function v that is 
equal to <f> on K, that vanishes outside the e-neighborhood of K, and has A k v G L 2 
for all k with \\A k v\\ 2 < 7c 1 c k C k t^ k e 2XT . 

On a first reading of the proof, we suggest thinking of the case where 0=1 and 
X has finite measure, so that Uj = 1 on Kj at each step. In this case the following 
heuristic may be helpful. 

Our goal is a function with two properties, the first of which is L 2 -smoothness, 
meaning that A k v £ L 2 for all k, and the second is the property of being = 1 
on K and = outside the e-neighborhood of K, which we call the characteristic 
property. Beginning with a function Ux which has the characteristic property but 
is not L 2 -smooth, we recursively apply a two step method. The first step smoothes 
iti by applying the heat operator for a small time i 2 to obtain v 2 = P t . 2 ui, which 
is smooth but does not have the characteristic property. The second step splits X 
into a neighborhood K% of K, and the complement Li of a larger neighborhood, as 
well as the region A2 between, and sets w 2 = 1 on K 2 , u 2 = on L 2 , and u 2 = u 2 
on A 2 . 

What we have gained in passing from m to it 2 is replacing the original abrupt 
drop of the characteristic function with the improved piece on A 2 , as illustrated on 
Figure [T] This argument is repeated inductively on nested annulus-like regions Aj , 
each time applying the heat kernel for a shorter time tj to get Vj = Pt Uj-i and 
then cutting Vj off to be constant outside Aj to get Uj . 

It is unsurprising that this process converges in L 2 . What is perhaps unex- 
pected is that A k v converges in L 2 for all k, and this is where the estimate (12. 3p 
is crucial. Essentially what is going on is that the "steepness" of the (j + l)-th 
interpolant depends both on the height it must interpolate and the "width scale" 
on which it interpolates. The "width scale" depends on tj+i through the norm 
II A k P tj+1 1| 2 ,2 < °ktj+i f° r the "steepness" measured by A fc , but the height to be 
interpolated depends instead on how much P tj changed the function during the 
smoothing step, which is controlled by £)(£,-, e2~ J ) because that is how much L 2 - 
norm can "leak" across the annulus Aj when we apply the heat flow. The series of 
terms A k (vj+i — Vj) is therefore essentially just tJ* 1 D(tj,e2~i), which converges 
by (EJl). 

Proof. Our assumptions imply that we can choose a sequence tj with tj = T 
and all tj < f such that (|2.3|) holds. 

Define neighborhoods Kj shrinking dyadically to K, neighborhoods Lj shrinking 
dyadically to the complement of the e-neighborhood of K, and annular regions Aj 
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Figure 1. Initial steps of the proof of Theorem 12.21 for a constant function. 

between them as follows: 

Kj = {x : dist(x,if) < e2- j }, 

Lj = {x : dist(x, K) > e(l - 2" J )}, 

Aj = {x : e2- j < dist(x, K) < e(l - 2~ j )} =X\ {K s U Lj). 
Using these regions, and writing Tj = ^ 2 *ti inductively define for j > 2 

ui = fXK x 
Vj = I', 1 

(P T J xeKj 

u j( x ) — \ v ji x ) x £ A/ 
[0 x e Lj 

where XKi is the characteristic function of the set K\. It will be notationally 
convenient to set t\ = and v\ = 0, so that the formula for Uj is valid for j = 1. 

From Lemma T2. 61 we see that Uj converges to a function v such that A fe w e L 2 
for all /c. However Uj and Uj have the same limit by Lemma 12.51 Since Uj = 
outside the e-neighborhood of K the same is true of v. Moreover, on K we know 
that 

u j = p Tjf -> 

and we conclude that the function v has the asserted properties. The estimate for 
A k v then follows immediately from (|2.6p and the estimate (|2.ip . because the latter 
gives 

||A fc i; 2 || 2 = |[P tl ui|| 2 < citf*||ui|| < ci^ fe . 

□ 

Lemma 2.3. For j > 2, 

\\(u j -v j ) X K j \\ 2 < £>(t i)e 2-')(l +11^-1 1| 3 ). 
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Proof. From the fact that Pt<fi = e~ xt cf> we have for j > 2: 
Pt.u.j ,-/V,/ 

On ifj the above is the difference between Vj and Uj. Since the support of the 
function to which we apply the heat operator is in Aj-% U £j— 1 and is therefore at 
least t2~^ distance from Aj we can use (12.2[) to see 



\ Uj -Vj\ 2 < D(tj,e2 3 ) 2 \\(Pr :i . 1 f(x)jx(x)A j . 1 uL j _ 1 + v j ^ 1 (x) x (x) Aj _ 1 \ 



\\( u j - v j)xk d \\ 2 < D(tj,e2 J )||(PT J _ 1 /(a;))x(a;)A J _ 1 uL J - 1 + v j -i{x)x{x)a^ 1 \\ 2 
<D(t 3 ,e2-i)(l + \\v J - 1 (x) X (x) A] _ 1 h) 
because Pt contracts L 2 norm and ||/||2 < 1- D 
Lemma 2.4. For j > 1 

j 

\\v j+ ih < ||«il|a < 1 + AY J D{t l ,,2- 1 ) < 3 

i=2 

Proof. Since P t is an L? contraction, the first inequality is immediate. Note that 
the second is true for j = 1 because ||ui||2 < H/H2 < 1- We also know from our 
assumption (|2.3p that £^ ,D(tj, e2~ J ) = C < 5. Induction and Lemma [231 then 
imply 

IMh = IbjXAj ||2 + \\UjXKih 

< \\VjXA, lb + \\VjXK, H2 + || («i - v 3 )xk 3 ||2 

< ||» J -|| a +X>(i i , e 2-^)(l+||t; J _ 1 || a ) 

< || Uj -_ 1 || 2 +£»(i J - I e2^)(l+lk--2||2) 

i-i i-2 

< 1 + 4 ^ D(ti, e2- 1 ) + D(t if e2-^)(2 + 4^ £2"*)) 

i=2 i=2 
J 

< 1 + 4^D(t 4 ,e2- ,; ) < 3 

i=2 

for j > 2. Of course when j — 2 the ||itj-2 1| term does not appear, and the estimate 
is trivial when j = 1 . □ 

Lemma 2.5. For j > 2 we have \\uj — Vj\\ 2 < 7D(tj,e2~i). 

Proof. Note first that Uj = on Lj, so 

\\( u 3 - v i)XL \\\ = \\v 3 XL \\\ = / \PtiUj-if < {D{t j ,c2- : >)f\\u j - 1 \\l 

because the support of is in Kj~i UAj—i, so is separated from Lj by a distance 
of at least e2~i , and therefore (|2.2[) is applicable. 
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Since Uj and Vj coincide on Aj we have 

ll u i _t, illa = U u i - v j)xk 3 \\ 2 + \\(uj -Vj)xl 3 \\ 2 

< D(i Jie 2^)(l + H^iX^-xlla) + (^(t J ,e2^))|| Uj _ 1 || 2 

< 7D(tj,e2- j ) 

by Lemma |2"31 and Lemma \2. 41 □ 
Lemma 2.6. For each non-negative integer k the sequence {A k vj} is I? -Cauchy. 
Proof. The heat operator commutes with A, so for j > 2 
\\A k (v J+1 - Vj )\\ 2 = \\A k (P t]+1 u 3 - Vj )\\ 2 

< \\A k P t]+1 (u 3 - Vj )\\ 2 + \\A k (P tj+lVj -vj)\\ 2 

(2.5) < Ic^D^^e) + At J+1 ||A fe ^|| 2 . 

where we used (|2.1[) . Lemma [2~5l and that ||Pf — /||2,2 < At for i < f by the spectral 
gap assumption. Then 

IIAS+1II2 < 70^0^,2-* e) + (l + A* i+ i)||ASlla 
and by induction 

||AS-+i|| 2 <7 C fcWv+i^ 2 ~' e ) II (l + AWi))||A fe t;2|| 2 

1=1 ^ m=/+l ' 

(2.6) <7c fe C fe e AT ||A fe t; 2 || 2 . 

because of (|2.3[) and the fact that J2j tj — T so Y\^°(l + \t m+ i) < e XT . Substituting 
into (|2.5p gives 

\\A k {v. ]+1 - v 3 )\\ 2 < TchtJ^D^tj^e) + 7c k C k \e XT \\ A k v 2 \\ 2 t 3+1 
which is summable by the assumption (|2.3p . □ 



The construction in Theorem 12.21 produces a function v satisfying A k v € L? for 
all k. The question of whether v is actually smooth (i.e. has A k v G C(X) for all 
k) is more difficult. Given that the only way we have to verify the assumptions 
for Theorem 12.21 in specific examples is to use an estimate of the form (|2.4p . and 
that in many instances the heat kernel is known to be continuous, the following is 
essentially as useful. 

Theorem 2.7. Suppose that the heat flow Pt is given by integration against a 
continuous kernel p(t, x, y) satisfying (|2.4|) . Let <fi £ A(X) and K be a fixed compact 
set. For any < e < eo there is a smooth function v that is equal to <j) on K and 
equal to zero outside the e-neighborhood of K . 



Theorem 12.71 is proved by applying Theorem 12.21 and the following lemma. 

Lemma 2.8. If Pt is given by integration against a continuous kernel p(t,x,y) 
satisfying ()2.4|) then A~ k maps L 2 (X,/i) into the continuous functions C(X) for 
all k > % — 1. In particular, if f is such that A k f £ L? for all k, then f is smooth. 
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Proof. The spectral representation of A implies 



A- k f(x) = 



t k p(t,x,y)dtf(y)dfi(y) 



t k p(t,x,y)dtf(y)dfi(y) 



however we may rewrite the second term using 



(t + l) k p(t + l,x,y)dt 



p(l,x,z) 

= (A y +I)- k p(l 1 x 1 y) 
By the inequalities of Holder and Minkowski 
A- k f(x) - A- k f(x') 



(t + l) k p(t + l,x,y)dtf{y)dn(y) 



p(l, x, z)p(t, z, y) dp,(z) dt 
(t + l) k p(t,z,y)dtdiJ,(z) 



< 



t k (p(t,x,y)-p(t, x',y))dt 



(A y + I)- k (p(l,x,y)-p(l,x',yf) 



< 



t k \\p(t,x,y) -p(t,x',y)\ 



L 2 (y) 



dt - 



L 2 (V) 



L 2 (y) 

+ as 



because (A y + I) k is a contraction of L 2 . Now \\p(t, x, y) — p(t, a; / ,y)|| i2 , . 

x — * x' , for any fixed t, because p(t, x, y) is continuous. This deals with the second 
term and shows that the integrand in the first term converges pointwise to zero. 
The latter is bounded by t k (\\p(t,x, -)||2 + \\p(t, x' , -)||2), which is an L 1 function 



when k > 





P(t,x,y) 



1, because 



p(t, x, y)p(t, x, y) dfi(y) = p(2t, x, x) < Cr a / p 



by (|2.4p . The result follows by dominated convergence. 



□ 



Corollary 2.9. If fi(X) < oo then for any compact set K and neighborhood U of 
K , there is a smooth non-negative function that is 1 on K and vanishes outside U . 

Proof. Apply the theorem to the constant function 1 and note that all of the Uj 
and Vj are non-negative, so the limit function is also. 

3. A SMOOTH BUMP AS A FIXED POINT OF AN OPERATOR 



□ 



To understand why it is sometimes possible to construct a smooth bump function 
on a self-similar set as a fixed point of an operator, we invite the reader to consider 
an elementary situation. Let / = [0, 1] be the unit interval in R. We may view 
/ as a p.c.f. self-similar set under the contractions fo — x/2 and fi — (x + l)/2. 
If /i is Lebesgue measure and £ is defined using a limit of a regular self-similar 
harmonic structures with resistance factors 1/2 then we obtain the usual Dirichlct 
energy and Laplacian, and the normal derivatives are the outward-directed one- 
sided derivatives at and 1 (see [T^l EE] for details). 
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The intuition for our construction is as follows. Consider a symmetric smooth 
bump function u on the interval / = [0, 1], for which u = 1 on [L, 1 — L] as shown 
in Figure O If we look at the graph of Aw = d 2 u/dx 2 we obtain something that 



FIGURE 2. The smooth bump function u. 



looks like a constant multiple of Figure 0J which appears as if it could be assembled 
from rescaled copies of u according to a rule like 



(3.1) 



H 



L 
2x-2 



if < x < I 
if j < x < L 
if L < x < 1 - L 
ifl-£<x<l-| 



1 



if 1 



< x < 1 



so that we might hope there is actually a smooth bump function u which has 
precisely this scaling behavior. If we let G denote the Green's operator for the 
operator A on J with Dirichlet boundary conditions, then this would be equivalent 
to asking that u be a fixed point of the operator 

G o $u(x) 



(3.2) 



^u(x) 



Go$u(l/2) 

It is a consequence of our general result Theorem 13.81 that the operator 'J in (|3.2 



1 



\ / 



L 



^-k, V 



Figure 3. The function Au = d 2 u/dx = $(u 



has a fixed point and that the fixed point is a smooth bump function. In fact more 
is true in the special case of /, where the fact that removing any interior point 
disconnects the set, along with the existence of an explicit formula for the Green's 
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function, allows us to prove that the fixed point has values in [0, 1] and is identically 
1 on [L, 1 — L]. For reasons of brevity we do not include the proof of this result; it 
is a simpler version of the proof of Theorem [ 



Proposition 3.1. If L is sufficiently small then the operator \1/ preserves the space 
of continuous functions on I that have values in [0, 1], vanish at and 1, and are 
identically 1 on [L, X — L]. Furthermore ^ is a contraction in the L°° norm on these 
functions, and its fixed point is a smooth function that vanishes to infinite order at 
and 1. 

Another example in which we can define operators $ and \& that are similar to 
(|3.ip and ()3.2p is the Sierpinski gasket SG with its standard harmonic structure 
and measure, where for sufficiently large I we can set 

f +1) {SG) 

(3.3) *u(x)=l-u(Ff 1 oFri(x)) 




if x e F> 

lix&FloF^SGlj^i 
otherwise 



and with p any vertex from V\ let 
(3.4) m u 



G o $u{x) 



G o $u(p) 

as illustrated in Figure [4] for the case 1 — 2. Again we omit the variant of the proof 




Figure 4. <Ht in the case I = 2. 



of Theorem 13.81 that establishes the following result 

Proposition 3.2. The operator of (|3.4[) is an L°° -contractive self-map of the 
set of functions that are continuous on SG, vanish at the boundary, are identically 
1 on SG \UiF,-, and satisfy \J u — 1 < i. The fixed point of ^ is a smooth bump 
function. 

The method described for I and SG rely heavily on the symmetry of these 
sets and on the assumption that they are endowed with the symmetrical harmonic 
structures and measure. This assumption is unavoidable if we want to use the same 
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operation <£> at all steps of the computation, as the natural linear combination of 
rescaled copies of the function will not otherwise have the desired properties, but 
it is very restrictive. Even some of the simplest of the nested fractals defined by 
Lindstr0m [17 have insufficient symmetry for a fixed $ to be used in the construc- 
tion of a smooth bump by this method. Nonetheless the method can be adapted to 
general p.c.f. self-similar fractals with regular harmonic structure and self-similar 
measure. 

Let X be p.c.f. self-similar with boundary Vq = dX, measure /x that is self- 
similar with scaling factors /ij and regular harmonic structure with factors Tj. We 
fix a scale l\ with size to be determined later, and label the boundary l\ cells by 
Yj = F^(X). Their union is Y = UYj. For any e > we will build a smooth 
function that satisfies \u — 1| < t on 1 \ 7 by a construction that inductively 
determines its Laplacian on the cells Yj , writing it as a fixed point of an operator 
^ on the following space of functions. 

Definition 3.3. Let C be the space of continuous functions u on X such that 
u(q) =0 for q 6 dX and \\u — l||i < \. Note that this space is non-empty and 
closed in the continuous functions with supremum norm. 

To define the operator "J we need a little more notation. Let S C VJ 1 consist of 
those points that lie in some Yj and in at least one other Zi-cell. If l\ is sufficiently 
large then no two of the Yj can intersect; we assume this and see that the connected 
components of X \ S are the cells Yj (less points of S) and the set X \ Y. Label 
those boundary points of the cell Yj at which Yj intersects another 1% cell by xt.j for 
i = 1, . . . , Ij. Fixing a second scale 1%, also with size to be determined, we associate 
to each Xij the unique (l\ + /2)-cell in Yj containing Xi j, calling it Z{ t j. We also 
set Zo t j = Fj 1+l2 (X), so it is the (l\ + ^-cell in Yj that contains qj € Vb, and 
define w^j to be the word such that Fij(X) — F Wi d (X) — Zij. Figure [5] illustrates 
our labelling conventions in the case X = SG, l\ = 2 and I2 = 1. We identify a 




9l 351,1 331,2 32 



Figure 5. Notation if X = SG, h =2 and h = \. 
particular function that is in C when l\ is large enough. Let / be the piecewise 
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harmonic function on X with values f{xi.j) = 1 for i = 1, . . . , Ij and j = 1, . . . , N 
but f(xo,j) — for all j. It is clear that / is continuous, identically 1 on X\Y and 
harmonic on each of the sets Yj. It fails to be harmonic only at the points £&y with 
i > 1, and we readily compute that the Laplacian of / is a measure supported at 
these points. In fact if S x denotes the Dirac mass at x then 

N Ij N Ij 

(3.5) A/ = a i>i 5 *i,i = ~ d nfj{xij)8 Xitj 

3=1 «=1 3 = 1 »=1 

with the second equality coming from the matching conditions for the Laplacian. 
Since fj o F* 1 is harmonic on X with boundary values at qj and 1 at all <?,•, i ^ j, 
it has normal derivatives bounded by a constant C(r) that may be chosen so as to 
depend only on the harmonic structure. Scaling of the normal derivatives under 
Fj 1 therefore implies \di,j\ — \®nfj(%i,j)\ < C(r)rJ Ll . 

The smooth bump function we seek will actually be a perturbation of /, con- 
structed by iteratively replacing the Dirac masses in (|3.5p by rescaled copies of the 
stage k bump, correcting for the boundary normal derivatives, and applying the 
Dirichlet Green's operator to obtain the stage k + 1 bump. We will see that each 
stage gains one order of smoothness, so the limiting function will be in dom(A°°). 
Our first step is to estimate the effect of a perturbation of the type described. 

Lemma 3.4. For each j = 1, . . . , N and i = 1, . . . , Ij let Vij be a finite, signed 
Borel measure with support in Zij . If we use the coefficients in (|3 . 5[) to define 

N Ij 
3 = 1 »=1 

and let u = G(v) be the result of applying the Dirichlet Green's operator, then 

(3.6) \u(x) | < C(r)N 2 sup ||^ j || for all x e X 

i,3 

where is the total variation ofvij. If in addition we have J = for all i 

and j then 

, N v 

(3.7) |u(i)| <C(r)N (sup \\ Vi j ||) ( ^ r k) for all x e X \Y 

hj \ fc=1 / 

Proof. Recall that G may be represented as integration against the continuous 
kernel —g(x,y), with sign chosen so g(x,y) > 0. The estimates we desire follow 
from (|1.4p and the fact that \aij\ < C(r)rJ l1 . The former ensures both that 
\9{ x iV)\ < C(r)rj on each Yj and that the oscillation of g(x,y) on Z^j is at most 
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C(r)r Wi j . We compute 

\u(x)\ = 



I g(x,y)dv(y) 

X 
N Ij 

<EEi fl ^ 

3=1 i=l 

N Ij 



\g{x,y)\dvij(y) 



7 = 1 i=l J Zi -i 



i(y)\ 



< 



C(r)N 2 sup \\vi j | 



which establishes the first inequality. To obtain the second we observe that Am = 
on X \ Y, so by the maximum principle we need only verify the inequality at the 
points Xij. Fix such a point j>, and use that each J dv^j = to subtract the 
appropriate constant from each integrand before estimating: 



Wxi',j')\ 



' g(xi>,j',y)dv{y) 

X 

N Ij 

{xi',j',y)dvi,j(y) 



E E 

3=1 i=l 
N Ij 

E E 

3 = 1 i=l 
N Ij 



>Zi 



< 



(g(xi',j',y) - 9{xi',j',x id fjdv itj (y) 

3 

J2J2 C ^ r J h / C{r)r Wi] d\ V% , 3 \ 
3=1 1=1 JZi -i 



N Ij 
3=1 i=l 



r k{i,j) ll^'jl 



because r 
fixed j , so we obtain 



= r}r£ for some k — k(i,j). Each k(i,j) occurs at most once for a 



,jv jv > , N 

(x VJ ,)\ < C(r)(sup H^ll) (EE^ 2 ) ^ <W(sup lk,||)(E 



□ 



There is an analogous but simpler estimate for the effect of introducing a mass 
supported on one of the small cells Zq j at the boundary. If vq j is a finite, signed 
Borel measure with support in Zqj we use (|1.4[) to see \g(x,y)\ < C(r)r' 1+ ' 2 on 
Zp.j and therefore 



(3.8) 



\G(v 0d )\ < 



' Zo.j 



\g(x,y)\d\v t) j(v)\<C(r)rf +l * HJ 



The ideas discussed so far allow us to generalize the definition of the operator 
^ in (|3.2j) and (|3.4[) . The idea is that to replace the Dirac mass terms on the 
boundary cells Z^j in (|3.5[) with normalized rescaled copies of u € C and apply the 



16 



ROGERS, STRICHARTZ, AND TEPLYAEV 



Green's operator to obtain a function that is near 1 on X \ Y. By adding some 
terms on the cells Zoj we can make the result have vanishing normal derivatives 
at the boundary without changing the value onX\y very much. In consequence 
we will obtain an operator that smooths u G C to be in C R dom(A) with vanishing 
normal derivatives and is near 1 on X \ Y . Iterating the operator will then produce 
a sequence of smoother and smoother bump functions. 
Let 

(3.9) u id {x) j, ■(./>//') '(«»^)) X x G Z itj 

I otherwise 

so that each Uij is continuous and has integral 1. Since u 6 C we also have that 
/K,| < (ju) \\u\\<3. 

Definition 3.5. The operator \& on C is tyu = G(v), where 

N Ij 

(3.10) v(x) X "<-^ n 

j=l i=0 

and G is the Dirichlet Green's operator. In this expression the coefficients for i > 1 
are given by bij — r/djj with aij as in (13. 5p . but the boj are yet to be determined. 



Note that \bi t j\ < C(r) when i > 1. It is immediate that G(v) = on dX and is 
continuous. Moreover AG(v) — v is a linear combination of continuous functions, 
so "fit G dom(A). The next lemma uses the Gauss-Green formula to reduce finding 
the correct &oj to a problem in linear algebra. 

Lemma 3.6. If l± and 1% are sufficiently large then there are values bo.j such that 
dn^u = on dX . The minimal sizes of 1% and li depend only on the harmonic 
structure of X and the number of vertices N in dX. 

Proof. Let hj be the function that is harmonic on X, equal to 1 at qj and at 
all other points of dX. Using AG(v) = v, G(v) e on dX, and the Gauss-Green 
formula 



E 

g k edx 



hjiy^d^y) = I fy(p)AG(«)(y)tfAi(») 

hj{qk)(d n G(v))(q k ) 

= {d n G(v))( qj ) 
from which d n G(v) = on dX is simply 

= / hj{y)v(y)dfj,(y) = ^b Vtj/ rj, h I hj(y)ui>j>(y)dfj,(y) 



for all j = 1, ... , N. Moving the terms depending on the fixed values bi/ji for i' > 1 
this may be reformulated as 

(3.11) 



J2 b oj' r j' h J hj{y)u ,f{y)dn(y) = - E E b i',f r j> h J 1 
i' ' -?'=! »'=i 
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which we recognize as a matrix equation ., M(u)j> t j 1 bo,j> = A(u)j with 



(3.12) 
(3.13) 



(M(ud/j,)) , 



3 ,3 
N !j' 



(A( U d^) 3 = - J2 E b 
j'=n'=i 



hj(y)uo,j'(y)d/j,(y) 

hj(y)ui> tj '(y)dfj,(y) 



r,3 jr j 



It is clear that we need to know M — M(udfi) is invertible, but rather than prove 
this directly we do so by proving a perturbation estimate similar to Lemma [3T4l that 
will be useful later. To this end consider replacing each of the measures uoj'd/j, 
in (|3.12p with a copy of a different probability measure da scaled and translated 
to give dai'j' supported on Zyji. We call the result M{da). The difference of 
these measures has mass zero, so we can compute an estimate involving the total 
variation of the measures 



M{udfi — da) j 



< 



hj(y)(u ,j'(y)dfx(y) - da ,f(y)) 

(hj(y) - hj(x ,f)) (u 0tf (y)dfi(y) - da ,j'(y)) 



(3.14) 



< C{r)r ll , +l2 \\ud^-da\ 



because hj is harmonic and therefore varies by at most C(r)r^ + ' 2 on each Zqji. In 
particular if the measures doy j> are Dirac masses at the points Xo,j* then M(da) 
is simply the identity, so (|3.14[) implies 

|(/-M)^| <C(r)^+ fa 

from which M is invertible when h+h is large, with M~ 1 \\ < C(N, r) J2 3 ; r / +2 - 
A similar perturbation argument can be made for A(dfi — da), where A(da) is 
obtained by replacing each Ui'j/dfj, by <iay j< in (I3.13| . Estimating the integral 
terms and using the bound < C(r) we obtain 



(3.15) 



\A{udii~da) j \ < C(N,r) (E r * 2 ) \\ ud V~M 



however this is not the most useful thing we can do here. Instead we recognize that 
the bounds \hj(y) — 1| < C(r)r 1 ^ on Yj and \hj(y)\ < on Yy for j' ^ j ensure 



A(ud[i)j 



E >■:>'' 3 

i' = l 



< C(N,r) 



so that combining this with our bound on / — M 1 we have 



(3.16) 



E^ 



< C(N, r)r 



If we examine the function / in (|3.5p it is clear that the normal derivative at 
each point xqj is Yli=i so our choice of bo j is a small perturbation of that 
which would be used to cancel the normal derivatives of /. We also remark that 
this shows all \b j\ < C(N, r). □ 
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If fa and I2 are large enough then the values bo.j from Lemma 13.61 may be used 
to complete Definition 13. 51 for the operator Some key properties of this operator 
are summarized in the following lemma. 

Lemma 3.7. If fa and fa are sufficiently large then ty(u) £ C PI dom(A) and 

(3.17) ll^l^ < d 

N 

(3.18) |*u(y) - l| < C 2 ^^ 2 forallyeX\Y 

where C\, C2 and the minimal sizes of h and I2 are constants depending only on 
the harmonic structure of X, the measure /i and the number of vertices N in dX . 

Proof. Since 

, N h N s 

*( m ) = G (J2J2 a ^ ui <j + b o,3 r j hu o,j ) 

S=l i=l j=l ' 

we obtain (j3~T7) from (|3l)j) and (pTg)) . and the fact that |6 ,j| < C(N,r) for all j. 
The estimate (|3 . 1 8|) is only a little more difficult. Using f(x) — 1 and (|3 . 5|) on the 

set X \ Y we see that 

|*u(x) - l| = \^u(x) - f(x)\ 

N Ij 



< 



, AT I j 
(■'{ ^ »/,, (''/..,<'/' - 

S=i i=i 



S=i i=l 

AT 



i=i 



(3.19) 



/ N \ ( N 

< C(N, r)8up||uijd/i- <^,J (J2 r k +C(N,r)(J2 
^ \=\ ' S=i 

<C(N,r)(j2^) 



rf \\u ,M\ 



and that for the aij terms is 



where the estimate for the bo ,j terms came from 

from p.7p because J Ui^dfi = 1 = J 8 Xi . and both are supported on Z tj3 . 
Finally we check that \\^u — l||i < \. Using the results we have so far 

||*u - l||i < C(N, r) (j2 r f ) <*M + j y (1 + Ci) dii 



(3.20) 



<C(N,r) (Y,r l f) + (1 + C 1 ) M (Y) 



so that we can be sure 6 C if both li and I2 are sufficiently large, because 
fi(Y) — > as li — » 00. It has already been observed that u € dom(A) and it e on 
dX, so the lemma is proven. □ 

Finally we come to the main result of this section. The following theorem im- 
plements the idea that motivated our definition of 'J, namely that ^> smoothes 
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functions in C and therefore its recursive application gives a bump function in 
dom(A°°). 

Theorem 3.8. Given e > there are l\ and I2 sufficiently large that ^ has a fixed 
point Uq in C with \u — 1| < e on X \ Y . The fixed point is a smooth bump function 
and every u € C has ||\& it — wo||oo — ► as k — > 00. 

Proof Let u,u £ C. We calculate — — G(v — 11), where v and v are as in 
(|3.10jl . Beginning with a variant of the computation ()3.19p we have 



\G(v-v)(x)\ < 



< 



(3.21) 



''•'f X X! •• "'^ 

S*=l i=l 

/ N h \ N 

S'=l i=l 7 i=l 

AT 

+ J2 r 7 h \ b o,j -bo,j\\G(uo,j){x)\ 
i=i 



A 



which suggests we will need to know estimates for both (ttjj — u^.j) and |6o,j — 6o,j | - 
Conveniently we can reduce the latter to the former using (|3.14p and (|3.15|) , because 
&0j and 6o,j are computed from equations of the form ^ ■ M(udfj,)j t ibojrJ 1 = 
A(udfi)i. We easily see that 

rj h (bo,j-bo,j) = ^M(ud^)~M(ud^-u^)^ + ^M(ud^) _1 M(ud^-ud^)(r~ /l 6 j)) 

however by (|3~T4")l we have both ||M(ttd/i-itd/z)|| < C(N,r)\\u - u\\ 1 J2 l r l l 1+l2 and 
that III-M-^ud/i)!! < c ( N i r )\\ u \\iY,i rl i 1+h , while (0113 gives us that P(ue^- 
< C(N, r) J^i r ' 2 — ""Hi- I n both cases we have used that the total variation 
of Uij — ttjj is bounded by |]itjj— ttjj ||i and that writing ux — J x u we can calculate 



\ui t j — Uij I = / u x 1 ?i(a;) — uxu(x) < u^} (l + lit^u^ 1 ) / l u — "I < — 



The conclusion is then that r rj ll \b j — b j | < C(N, r)\\u — u\\i J2i r ; 2 ■ Substituting 
this into (|3~2Tj) and using IpTTj) and (f3T8]> we find that on X \ Y 

|G(v-v)(a:)| < C*(7V, r) ( ^ r' 2 J sup|[ti i(J -d/x - u^d^H + c Wr)(^\ 

^ 4 — 1 



r{ u-u 1 



because the total variation — Ujj<i/i|| was already computed to be at most 

8||u— On the rest of X we must use (|3.6p instead of (|3.7p . The weaker estimate 
is easily computed to be 



(3.22) 



\G(v-v)(x)\ < C(JV,r)||« — 
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From our estimates on G(v — v) = — f m we see that ^ is a contraction on 
L 1 if l\ and I2 are sufficiently large, because 

Htftt-*^ < fj,(X\Y)C(N,r)(y2r l ^\\u-u\\ 1 + (i(Y)C(N,r)\\u - u\\i 

i 

< C(N, r) + r \ ? ) h-^\\i 

3 

It follows readily that iff has a unique fixed point in C and iff k u converges to this 
fixed point in L 1 . From (|3 . 22[) this convergence is uniform, and we notice that the 
correct choice of I2 provides |u — 1 1 = \iffuo — 1| < e on X \ Y by (I3.I8|) . 

It remains only to see that uq is a smooth bump function on X. Inductively sup- 
pose f! k u e dom(A fc ) and both A^ k u = on dX for < j < k and d n A^ k u = 
on dX for < j < k — 1. This is certainly true for k = 0. By construction, A^f +1 u 
is a linear combination of rescaled copies of iff k u that have been extended by zero 
as in (|3.9p . Each of these functions is in dom(A fc ) by the matching conditions for 
the Laplacian, so we conclude that \E' fc+1 u € dom(A fc+1 ). It is immediate that 
A^ fc+1 u = on dX for 1 < j < k + 1 and d n A^ k+1 u = on dX for 1 < j < k. 
By Lemma [3~71 we know also that fy k+1 u and d n *5> k+1 u vanish on dX, which closes 
the induction and establishes that uo G dom(A°°) and vanishes to infinite order on 
dX. □ 

4. A BOREL THEOREM ON P.C.F. FRACTALS 

The classical Borel theorem tells us that given any neighborhood of Xo £ R and 
any prescribed sequence of values for u and its derivatives at Xq , we may construct 
a smooth function u with support in the neighborhood and the given sequence 
of derivatives at Xq. Using the smooth bump functions we have constructed, we 
now show that the same result holds at junction points of certain p.c.f. fractals. 
In what follows X is p.c.f. and self-similar under {Fj}^L 1 and the measure fj, is 

self-similar with factors < fij < 1, Yli = 1' so tnat ^{^w{X)) — JXjLi ^wj 
when w is the word w\ . . . w m . The Dirichlet form is that associated to a regular 
self-similar harmonic structure with resistance renormalization factors < rj < 
I for j = 1,...,N. Our arguments depend on the existence of smooth bumps 
as previously constructed. The crucial step is the existence of smooth functions 
with finitely many prescribed normal derivative values, which is established in the 
following lemma. 

Lemma 4.1. Given a boundary point q £ Vq there are smooth functions /; such 
that 

A k fi(p) = for allk>0 and all p £ V 
d n A k fi(p) = for allk>0 and all p£V \ {q} 
d n A k f l (q) = S lk 

Proof. We begin with the case I — 0. If U is the smooth bump function on X 
produced in Theorem 13.81 we localize it near the boundary points of X at a scale 
to to be determined later. Define 

u . = l UoF f m ™Fp{X) 
1 otherwise 
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and observe from the matching conditions for the Laplacian that each Uj is smooth. 
Now apply the Dirichlet Green's operator G to these functions and form the linear 
combination 

JV 

with coefficients to be chosen. It is clear from the properties of U that A k f = 
on Vq for all k > and that d n A k f — on Vq if k > 1. Moreover the Gauss-Green 
formula yields values of the normal derivatives at the points q t G Vq 



N 

E« 



h„U, 



x 



where hi is the harmonic function on X with boundary values hi(qj) — u tJ . 
then suffices that we can invert the matrix with entries Aij = J hiU : 



Sij. In 

order that there be coefficients dj such that / has the properties asserted for /o it 

j. We use the 
fact that 

\hi\<rf on Ff(X) for j^i 

|Ai-l|<rr on F™{X) 

which follow from the estimates on the oscillation of a harmonic function that were 
mentioned at the end of the introduction. Using this we calculate 



(4.1) 



hiUi 



U 



\U\ fovj^i 



X 



X 



(hi -l)Ui 



\U\ 



x 



Let D be the diagonal matrix with entries Da = /i™ J x U . Then we readily compute 

(AD^ 1 )^ = fij m (f ll) 1 Aij is close to the identity if m is large. Indeed, by (|4.ip 

we have \(I - AD- 1 )^ < Cp m with p = maxir, and C = (J Uy\f \U\), so that 
AD^ 1 is invertible provided m is sufficiently large. 

We proceed by induction on I, with an almost unchanged argument. Suppose 
the functions /; for I < L — 1 have been constructed as linear combinations of the 
form 

/ N 

(4.2) f l = Y J T. a ^G {n+1 \U ] ) 1<L-1 

n=l j=l 

and consider the function 

N 
3 = 1 

so that A k f = on V for all k and d n A k f - on V for k > L + 1. When k = L 
we have 



A 



d n A L f( qi 



a i Ai - 



where Aij is as before, so we may select a, to obtain d n A L f(q) = 1 and d n A L f(p) = 
at other points p G Vq. Subtracting an appropriate linear combination of the // 
for I < L — 1 we obtain the desired /l in the form of (|4.2[) . □ 
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With this in hand it is simple to deal with finitely many values of the Laplacian 
at a boundary point q GVo- 

Lemma 4.2. Given a boundary point q £ Vq there are smooth functions gi such 
that 

d n A k gi(p) — for all k > and all p e V 
A k gi(p) = for all k > and all p £ V \ {q} 
& k gi(q) = Si k 

Proof. Let h be the harmonic function which is 1 at q and at all other points of 
Vo- Clearly A k h = for all k > 1 and therefore also d n A k h(p) — for all k > 1 and 
p G Vo- For each p £ Vq let /o lP be the function constructed in Lemma l4.1l with non- 
vanishing normal derivative at p. It is clear that go — h— XLgvb (pnh{p))fo,p has the 
desired properties, so we have found the first of our functions. To obtain the others 
we simply apply the Dirichlet Green's operator G. Notice that A k G l h(p) — SkjS p ^ q 
for all k and all p £ Vq, and also d n A k G l h(p) — for all k > I + 1. To obtain it 
remains only to subtract off all normal derivatives that occur for < k < I using 
the functions from Lemma l4.ll □ 

The proof of a Borel-type theorem from the above lemmas is standard. All 
that is needed is information about how scaling the support of a function changes 
its Laplacian and normal derivatives. Recall that for a Dirichlet form associated 
to a regular self-similar resistance, both the Laplacian and the normal derivative 
may be obtained as renormalized limits of corresponding quantities defined on the 
approximating graphs (Section 3.7 of [14j). In particular, pre-composition with 
the map F~ l rescales the fc-th power of the Laplacian by {^iri)~ k and its normal 
derivative by n~ k . For this reason, if q = qi is the boundary point of interest 
we define 



fl 



l.m 



9l, 



> ^l^d+l) fioF (-rn) QnFr 

otherwise 
otherwise 



so that we have for all k 

A%n(3)=0 A k g Lm (q) = 5i k 

(4.3) d n A k f Lm (q) = S lk d n A k 9l , m (q) = 

but the L°° norms of the lower order derivatives have decreased and those of the 
higher order derivatives have increased. 

(4.4) ||A fe /, m ||oc = ^ k) rT (l+1 ~ k) \\^fi\\oo 

(4.5) l|A fc fl!, m ||oo = (^r (Z - fe) ||A fc 5i ||oo 

With this in hand we can easily prove our version of the Borel theorem. 

Theorem 4.3. Let q £ Vq be fixed, and be an open neighborhood of q. Given a 
jet p = (pq , pi , . . . ) of values for powers of the Laplacian and a — (<7o of 
values for their normal derivatives, there is a smooth function f with support in O 
and both A k f(q) = p k and d n A k f(q) = a k for all k. 
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Proof. We give the usual proof that it is possible to define / by the series 
( 4 -6) / = ( Pl9l,m t + o-ifi <n , 



for an appropriate choice of mi and n;. 

Let mo = no be sufficiently large that F^^X) C 57. For each I choose mi > mo 
so large that 

Wp^gi^W^ < 2 k - 1 - 1 for < k < I - 1 

using the scaling estimate (|4.5|) . Similarly use the scaling relation (|4.4ll to take 
ni > mo such that 

||^A fc /i,n,|| 00 <2 k - 1 - 1 forO<fc<?-l 
Then for fixed k we have 



oo 



l=k+l 



< ^2 k ~ l < 1 



fc+i 



so that all powers of the Laplacian applied to (|4.6[) produce L°° convergent series. 
It follows that / as defined in (|4. 6|) is smooth and has support in f2. By (|4.3|) it 
has the desired jet, so the result follows. □ 

We remark that for any e > we could replace the bounds 2 fc ~'~ 1 in the proof 
with e2 k ~ l ~ 1 . It follows that we can define / by (14. 6|) and have the estimate 



(4-7) ||A fc /L < C(k,n)J2(\Pi\ + Wi\) +e 

where C(k, f2) does not depend on the jet we prescribe. 

It is also possible to estimate the effect of the size of the support F/™ (X) on a 
finite collection of jet terms. Take mi = mo for < j < L, ni — mo for < j < L— 1 
and thereafter choose both so large that the estimates e2 k ~ l ~ 1 hold. From (|4.4p 
and (|4.5p we deduce for < k < L, 

(4.8) iiA fe /iu < c(fc)(r^r mofe (j2( r ^r al \pi\ + E r r o(i+1 VrVd) + ^ 

Similar estimates hold for other norms and seminorms, provided only that their 
scaling is understood. For example, if we consider the seminorm £{A k f)^ 2 , then 
the scaling in (|4.4[) and (|4.5|) is reduced by r~ m ^ 2 , so the above construction would 



give 

£(A k f) x/2 < c(fc)rr m ^ fc+1 /2) M -™o^^ (rt/ii) ™o^ i | + 



The result of the theorem may be transferred to any junction point F w (q) and 
cell F w (X) in X, simply by modifying the desired jet to account for the effect of 
composition with F w , solving for / on X , and defining the new function to be 
/ o F^ 1 on the cell. We record a version of this that will be useful later; note that 
in the following we use the notation d n for the normal derivative with respect to 
the cell F W (X). 
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Corollary 4.4. Let F w (q) be a junction point in X. Given a jet (po, pi, ■ ■ 
(ci, 02, . . . ) there is a smooth function f on F W (X) that has A k f{p) — d n A k f{p) = 
at all points p s dF w (X) such that p ^ F w (q), and satisfies A k f(F w (q)) = pk 
and d n A k f(F w (q)) = ak for all k. 

5. Additive Partitions of Functions 

The results of [4] show that multiplication is not generally a good operation 
on functions in dom(A). In particular, for X a p.c.f. fractal with self-similar 
measure and regular self-similar harmonic structure it is generically the case that 
if u € dom(A) then u 2 £ dom(A). In such a situation there is no hope of using 
a smooth partition of unity to localize problems in the classical manner. Instead 
we provide a simple method for making a smooth decomposition of / G dom(A°°) 
using Theorem 14.31 Throughout this section we make the same assumptions on X 
as were made in Section SJ 

Theorem 5.1. Let U Q f2 Q be an open cover of X and f £ dom(A°°). There is a 
decomposition f = Y]u—i fk in which each k has a corresponding ak such that fk 
is smooth on X and supported in Q ak . 

Proof. Compactness of X allows us to reduce to the case of a finite cover uf fl ak 
for which there is no sub-collection that covers X. We write f2fc for £l ak , and 
construct the functions fk inductively. At the fc-th stage we suppose there are 
functions f±, .. . , fk-i with the properties asserted in the lemma, and that function 
gk-i = f — Y]i—i fi is smooth on X and vanishes identically on a neighborhood 
Tlk~i of X \ (u?_ k Qj\ ■ In the base case k = 1 this assumption is trivial, and it is 
clear that the theorem follows immediately from the induction. We have therefore 
reduced to the case where our cover consists of the two sets fife and f2/- = U^L A , +1 f2 J -, 
because the induction is complete once we have fk as in the lemma such that gk is 
identically zero on a neighborhood Ilk of X \ Sl^. 

For a scale m and x € X, define the m-scale open neighborhood of x to be 
the interior of the unique m-cell containing x if x V m , and to be the union 
{x}U (u w F w (X)\dF w (X)) if x = F w (qi) is a junction point. By Section 1.3 of [H] . 
the m-scale open neighborhoods form a fundamental system of open neighborhoods 
of x. At each x in fife there is a largest m such that the m-scale neighborhood of 
x is contained in f^, . The collection of all such largest neighborhoods of points of 
ilk is an open cover of the compact set Sppt((/fc_i) \ (lk- We use to denote the 
union over a finite subcover. 

Clearly has finitely many boundary points. Let those boundary points that 
are also in &k be Xi, . . . ,xj, and take at each a finite collection of cells {Cij}^ 
having Xj in their boundary. We require that U ^U^CVj^ contains a neighbor- 
hood of Xj, that all of the Cij lie entirely within Qk an d none intersect Ak, and that 
Cjj H Cy ji is empty unless j — j' , in which case it contains only xj. On each cell 
we apply Corollary 14.41 to find functions hij that match gk—i(xj) and all powers 
of its Laplacian at Xj, and such that the sum J^. hi j has normal derivatives that 
cancel d n A n gk-i(xj) at Xj for all n. Thus matches gk-i in the sense of 

the matching condition and vanishes to infinite order at the other boundary points 
of U^Cij. The matching condition implies that fk — <?fc-i|A fc + J2jJ2i^i,j * s 
smooth. It is clearly supported on and equal to / on the closure A& of Ak, 
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so <7fc = gk-i — fk is zero on a neighborhood Hk of X \ Qk, which completes the 
induction and the proof. □ 

Remark 5.2. Theorem 15.11 is stronger than the partitioning results we obtained 
directly from the heat kernel smoothing procedure of Section [21 but only applies to 
the more limited case of p.c.f. self-similar sets with regular Dirichlet form. We note, 
however, that it provides a stronger version of Theorem 13.81 because by smoothly 
cutting off the constant function 1 we obtain a smooth bump that is identically 1 
on a compact K and outside an e-neighborhood of K . In contrast to Corollarv l2.9l 
it does not imply that there is a positive smooth bump of this type. 
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